Field induced transition of the S = 1 antiferromagnetic chain with anisotropy 
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The ground state magnetization process of the S — 1 antiferromagnetic chain with the easy-axis 
single-ion anisotropy described by negative D is investigated. It is numerically found that a phase 
transition between two different gapless phases occurs at an intermediate magnetic field between 
the starting and saturation points of the magnetization for —1.49 < D < —0.35. The transition 
is similar to the spin flopping, but it is second-order and not accompanied with any significant 
anomalous behaviors in the magnetization curve. We also present the phase diagrams in the m-D 
and H-D planes which reveal a possible re-entrant transition. 

PACS Numbers: 75.10.Jm, 75.40. Cx, 75.45. +j 



I. INTRODUCTION 

The one-dimensional (ID) antiferromagnets show var- 
ious strange phenomena due to quantum effects. One 
of interesting examples appears in the low-temperature 
magnetization process of anisotropic antiferromagnets 
with an easy axis for the Neel order. In higher dimen- 
sion including 2D the external magnetic field H along 
the easy axis gives rise to a first-order phase transition 
called spin flopping which results in a jump in the mag- 
netization m of the S — \ Heisenberg antiferromagnct 
with the Ising-like anisotropic coupling. JlJ In ID, how- 
ever, the transition was shown j^] to be a second-order 
one described by the critical behavior m ~ \H — H c \ x / 2 , 
where H c is the critical field. It implies that the discon- 
tinuity of the magnetization curve vanishes due to the 
strong quantum effect in ID. 

The S = 1 systems can have another-type easy-axis 
anisotropy described by Dj2j(Sj) 2 (D < 0) which 
would also yield the spin flopping in higher dimension. 
The ID S — 1 antiferromagnet, ,even in the isotropic 
case, has a quite different feature from S = h, which is 
characterized by a gap in the low-lying excitation spec- 
trum, called Haldane gap (HG). j| Thus the magneti- 
zation process of the S = 1 chain with the easy-axis 
anisotropy could include some interesting problems. In 
contrast to the easy-plane anisotropy (D > 0), the nega- 
tive D case has been investigated in few works, partially 
because a noted quasi-lD S = 1 antiferromagnet with 
negative D CSN1CI3 was revealed to have the Neel order 
at low temperatures, which implies it doesn't have the 
Haldane gap. However, it could be a good material to 
test the spin flopping problem. 

The ground state (GS) magnetization curve of the ID 
5 = 1 isotropic antiferromagnet has already been given 
using the finite cluster calculation and size scaling |f|] 
based on the conformal field theory(CFT). || It has two 
critical fields H c \ and H C 2 which denote the starting and 
saturation points the magnetization, respectively. H c \ 
corresponds to the amplitude of HG, and GS has no 
magnetization and the exponentially decaying spin cor- 



relation for H < H c i. The size scaling analysis also sug- 
gested that the transitions at H cl and H c2 are second- 
order ones with the critical behaviors m ~ (H — Hd) 1 ' 2 
and 1 — m ~ {H C 2 — H) 1 / 2 , respectively. The further 
analysis || on the excitation spectrum revealed that for 
H c i < H < H C 2 the low-energy collective mode of the 
quasiparticles created by SjS is well described by the 
interacting fermion system called Luttinger liquid (LL) 
which is characterized by the gapless charge and cur- 
rent excitations. In terms of the original spin system it 
implies that the excitation changing the uniform magne- 
tization and the soft mode of the incommensurate spin 
density wave with the wave vector 2kp = 2nm are gapless 
for < m < 1. 

In the nonmagnetic GS of the ID S = 1 antiferro- 
magnet the negative D enhances the antiferromagnetic 
spin correlation along the easy axis and the system has 
the Neel order for sufficiently large D. The phenomeno- 
logical renormalization group (PRG) analysis || up to 
the system size L=20 gave the Ising-like critical point 
D c i = —0.35 ±0.03 which separates HG and Neel phases. 
For HG phase the magnetization process is expected to 
be qualitatively equivalent to the isotropic case. Note 
that even if HG vanishes due to the negative large D, 
H c i is still finite because the magnetic branches with 
^2S Z = ±1 of the first excited triplet state go up in 
the energy spectrum due to the negative D, in contrast 
to the positive D. Since large negative D excludes the 
state with = at every site, each spin behaves like 
S = 2- I n f ne limit D — > — 00 the system is equiva- 
lent to the antiferromagnetic Ising model which shows 
only a trivial magnetization process described by only 
one step at H c \ = H C 2. The second-order perturbation of 
the exchange interaction, however, leads to the effective 
exchange process between | — 1,1) and [ 1, — 1) at adja- 
cent sites. Thus the S = 1 chain with large but finite 
negative D is expected to be equivalent to the S = ^ 
chain with Ising-like anisotropic couping. The equiva- 
lence seems to suggest that the S = 1 chain with large 
negative D shows quite similar magnetization process to 
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the isotropic case; the two second-order phase transitions 
with the same critical behaviors as the isotropic case oc- 
cur at H cl , which corresponds to the Ising gap, and H c2 . 
The low-lying excitation of the system is also predicted 
to be described by LL for H c i < H < H C 2. The charac- 
ters of the gapless excitations, however, arc different from 
the isotropic case. Since the state with 5| — leads to a 
large energy loss proportional to D, the quasiparticle of 
the gapless magnetic excitation should be described by 
{Sj^) 2 instead of Sj. In addition the momentum of the 
soft mode is known to be 2kp = (1 — to)tt, while it is 
2kp — 2-Km for the isotropic system. In terms of the in- 
teracting fermion systems, the difference in the soft mode 
is attributed to the change of the band picture which 
results from the pair creation due to the large on-site 
attraction made by D. Thus we distinguish LL phase 
that appears in the large negative D region, denoted as 
LL2, from another LL phase including the isotropic case, 
denoted as LL1. 

In this paper we investigate the phase boundary be- 
tween these two LL phases using the finite cluster calcu- 
lation up to L — 20 and some size scaling techniques. In 
addition we consider the possibility of the field induced 
phase transition, based on some phase diagrams, and the 
relation to the spin flopping. 

II. MODEL AND NUMERICAL CALCULATIONS 



the D-J z plane, where J z is the z-component of the ex- 
change coupling, obtained by a bosonization technique. 
It suggested that the phase boundary was Ising-like. 
Consulting the result from the bosonization, to clarify the 
features of LL1 and LL2 in GS, we consider the following 
three spin correlation functions: 

C 1 (r)=(S +5-)~(-irr-" 1 , 
C 2 (r)H(S+) 2 (S-) 2 )~r-^, (2) 
C z (r)=(S z S^) ~cos(2k F )r-^, 

where we take only the most dominant term for each. 
Ci(r) and C2(r) are associated with the gapless single- 
particle (SM = ±1) and two-particle (5M — ±2) exci- 
tations. C z (r) is related to the gapless soft mode with 
2k p. It is expected that both Ci(r) and C2(r) obey the 
power law shown in (j|) in LL1, while only C\(r) decays 
exponentially in LL2 because the pair excitation is gap- 
less, but the excitation with 5M = ±1 is massive in LL2. 
It is also predicted that the momentum of the soft mode 
corresponds to 2kp i = 2~Km for LL1, 2kp2 = (1 — m)7r 
for LL2. Therefore, if we consider the single-particle ex- 
citation gap Ai and the soft mode gap with the momen- 
tum 2kpi A2fe F2 , the two gaps should cross at the phase 
boundary between LL1 and LL2. Because Ai = and 
A 2feF2 ^ in LL1, while Ai ^ and A 2fcF2 = in LL2 
in the thermodynamic limit. 



We start with the ID S = 1 antiferromagnetic Heisen- 
berg Hamiltonian with the single-ion anisotropy in a 
magnetic field 

H = Ho + Hz, 

H =£ s ,- s ,+i + ^E(^) 2 ' (i) 

j j 

Hz= —H^y^s*, 

j 

under the periodic boundary condition. We consider only 
the easy-axis anisotropy (D < 0) and the Zeeman term 
Hz implies the applied magnetic field along the axis. For 
L-site systems, the lowest energy of Ho in the subspace 
where the eigenvalue of J^j &j i s M (the macroscopic 
magnetization is m = M/L) and the momentum is k, is 
denoted as Ek(L,M). In addition we define E(L,M) as 
the lowest one among Ek(L, M)'s. Using Lanczos algo- 
rithm, we calculated E k (L,M) (M = 0,1,2, ••-,L) for 
even-site systems up to L — 20. 



III. SPIN CORRELATION FUNCTIONS AND 
SOFT MODES 

The two phases LL1 and LL2 are similar to the two 
planar phases in the nonmagnetic GS phase diagram in 
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FIG. 1. Scaled gaps LAi (x), LA^ /2 (o) and LA 2 (A) are 
plotted vs D for L = 20, 16 and 12 with fixed m = 1/2. LA 2 
does not depend on L, which means it is always gapless in the 
thermodynamic limit. The crossing point of LAi and 
is almost independent of L and gives a good estimation of the 
critical point as D c — —1.44 ± 0.01. 



IV. PHASE BOUNDARY AND CRITICAL 
INDICES 

In order to investigate the phase boundary between 
LL1 and LL2 in GS with m ^ 0, we vary the parameter 
D(< 0) with fixed m — 4£. PRG is a good approach to 
study on phase transitions, particularly with the Ising- 
like universality. If we assume the critical behavior of the 
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gap A ~ \D — D c \ v , the L-dependent fixed point D C (L) 
can be determined by the PRG equation 

L L+SL A(L + 5L,D r ) = LA(L,D), (3) 

and the L-dependent exponent v{L) is also estimated by 



v(L) = In [ 



L + 6L- 



/In 



- (L + SL)A'(D C (L),L + 6L) - 
■ LA'(D C (L),L) - 



(4) 



where A'(L>, L) is the derivative of A(D, L) with respect 
to D. The extrapolation of D C (L) and v{V) to L — > oo 
give D c and ^ in the thermodynamic limit. 

Based on the above argument, the two gaps Ai and 
A2fe F2 are taken for the order parameters. The two gaps 
of the finite size systems are calculated by the forms 

Ai = E n (L, M + 1) + E«(L, M - 1) - 2E Q (L, M), 
A 2kF2 = E 2kF2 (L,M) - E (L,M), (5) 

where M — mL, and we restrict M to an even number 
because odd M states cannot be the ground state in LL2. 
The form of Ai is obtained by the sum of the two gaps 
with 5M = +1 and —1 of the total Hamiltonian H, to 
cancel out the Zeeman energies in both gaps. 

For to = | where 2kp = -| for LL2, the scaled gaps 
LAi and LA n / 2 are plotted versus D for L =20, 16 and 
12 in Fig. [I]. The two-particle excitation gap 

A 2 = E Q (L, M + 2) + E (L, M - 2) - 2L (L, M), (6) 

is also plotted as LA 2 in Fig. [l[ LA 2 is almost inde- 
pendent of L in the whole region, which means the two- 
particle excitation is always gapless because the size de- 
pendence of the gap is A ~ at a gapless point, as CFT 
predicted. On the other hand, the behaviors of Ai and 
Att/2 obviously suggest the existence of a critical point 
D c between LL1 and LL2. They imply that Ai (A„./ 2 ) is 
gapless (gapped) in LL1, while gapped (gapless) in LL2. 
In Fig. [ythe L-dependent fixed points of Ai and A^m, 
in addition the crossing point of the two gaps with fixed 
L, look almost independent of L. The extrapolation of 
the both L-dependent fixed points for SL = 4 using the 
standard least-aquare fitting of the quadratic function of 
gives a common value D c = —1.45, which shows a 
good agreement with the extrapolated result of the cross- 
ing point D c = —1.44, as shown in Fig. ||(a). Thus we 
take the result of the crossing point of Ai and A^/ 2 f° r 
the best estimation as D c = —1.44 ± 0.01. The same ex- 
trapolation of the exponent v gives the conclusion v = 1 
for both gaps with the error less than a few percent, as 
shown in Fig. |^(b). It is consistent with the Ising uni- 
versality. 
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FIG. 2. (a)Some estimations of the critical point D c for 
m = 1/2: the L-dependent fixed points D C (L) for Ai (solid 
circle) and A w /2 (solid diamond), and the crossing point of 
the two gaps (solid triangle) plotted versus 1/(L + 2). The 
extrapolation by fitting of the quadratic function of l/(L + 2) 
gives D c = —1.45, —1.45 and —1.44 ± 0.01, respectively. 
(b)L-dependent exponents v(IS) for Ai (solid square) and 
A^/2 (solid triangle down) plotted versus 1/(L + 2). The 
same extrapolation gives v — 1 for both gaps with the error 
less than a few percent. They are consistent with the Ising 
like universality. 

For further test of the equivalency of the present phase 
boundary between LL1 and LL2 to the one between the 
two planar phases in the nonmagnetic ground state phase 
diagram discussed above, we estimate the exponents rji, 
r\ 2 and r\ z for 2kp 2 in the form (Q), using the L depen- 
dence of the gap associated with each spin correlation 
function like 



A ~ nv s rj/L, 



(7) 



where v s is the sound velocity which is obtained from the 
form 

^-(E 27r/L (L,M) - E (L,M))+O(^), (8) 



as CFT predicted. For m — i the estimated expo- 
nents rji, rj 2 and r\ Zl using the extrapolation of the 
size-dependent results from the form (0) by fitting the 
quadratic function of are shown in Fig. |3| Except 
for the rapidly diverging behaviors of r\\ and r\ z near the 
critical point D c = —1.44 where the size correction is 
too large to estimate them, all the exponents seem to ap- 
proach the expected values 771 = j, r\ 2 = 1 and r) z — 1, 
with D approaching D c . In addition the values of 4771/772 
and rj z rj 2 are plotted in Fig. |j, to check the universal 
relations 772 = 4?7i for LL1 and r\ z r\ 2 = 1 for LL2. The 
relations are revealed to be satisfied at least as far from 
the critical point as \D — D c \ > 0.3. They are all con- 
sistent with the results of the bosonization for the two 
phases in the nonmagnetic ground state. Our estimation 
of the central charge c of CFT using the form of the size 
dependence of the ground state energy 



1 7T 1 

-E (L,M) - e(m) - qCv s j^, 



(9) 
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results in c = 1 at least sufficiently far from D c , as ex- 
pected for both LL1 and LL2, while c — 2 just at D c . 
It implies that two Ising-like criticalities with c = | for 
each described by Ai and A2fe F2 appear at D c on the 
background with c = 1. We also got the same result of 
every critical exponent for m = 4 and |. 



Dc 




-2.0 -1.5 -1.0 -0.5 0.0 

D 



FIG. 3. Estimated exponents 771 (dotted line), 772 (solid 
line) and r\ z (dashed line) are plotted versus D for m = 1/2. 
The values of 4771/772 (long dashed line) and 772772 (dot-dashed 
line) are also plotted to test the universal relations 772 = 477 
for LL1 and 772772 = 1 for LL2. 



V. PHASE DIAGRAMS 
A. m-D phase diagram 

The crossing point of the two gaps Ai and A2k F2 > f° r 
various magnetizations that can be obtained by m = j 
using even M for L=16, 18 and 20, are plotted in the 
m-D plane in Fig. ^. Since they give a common curve 
almost independent of L, we don't consider any size cor- 
rections except for the limits to — > and 1. The ex- 
trapolation of the crossing point for various L with fixed 
M = 2 leads to the phase boundary in the limit m = 0+ 
as D c i = —0.35 ± 0.01, which gives a good coincidence 
with the boundary between HG and Neel phases in the 
nonmagnetic GS derived from PRG. A similar estimation 
yields D c2 = -1.32 ± 0.01 in the limit m = 1-. Fitting 
a suitable polynomial to D c for L = 20 and 18, and the 
extrapolated D c \ and D C 2, the conclusive phase bound- 
ary is obtained as a solid line in Fig. |]. It suggests that 
the phase transition from LL2 to LL1 occurs in the mag- 
netization process for D c % < D < D c \, where D c % is the 
minimum of D c , as shown in Fig. ^, and D c3 = —1.49 
is obtained by the fitting curve. In addition the phase 
diagram shows the existence of a re-entrant transition 
returning to LL2 for D c3 < D < D C 2- 
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FIG. 4. m-D phase diagram. The boundary between 
LL1 and LL2 is determined by the polynomial curve fitting 
(a solid line) to the crossing points of LAi and LA2k F2 f° r 
various m's with L = 20 and 18, which are almost indepen- 
dent of L. The extrapolations of D c for M = 2andL — 2 gave 
D c i = -0.35 ± 0.01 and D c2 = -1.32 ± 0.01, respectively. 
The above curve fitting yielded D C 3 = —1.49. In addition the 
classical spin flopping line is shown as a dot-dashed line. 

B. H-D phase diagram 

Since the excitation with SM — ±1 is gapless for < 
rn < 1 in LL1, CFT predicts the asymptotic form 

± E(L, M ± 1) T E(L, M) - H ± Ay + o(-) (10) 

where A = nv s rji. Thus the average of the two forms 
(pp|), canceling the dominant size corrections, gives a 
good estimation of H for to = ^ as 

H$ = [E(L, M + 1) - E(L, M - l)]/2 = H + o(i). (11) 

To get the phase boundary between LL1 and LL2 in H-D 
plane, we calculate at D c for given to with L = 20 

and 18, and plot (D c ,h9j) f° r various to's. The curve 
obtained by the polynomial fitting to the plot, is shown in 
Fig. ||, where H c3 is the lower boundary and H ci is the 
upper. H c i is estimated by the Shanks transformation 
fUJ applied to E(L, 1)-E(L, 0) for L = 12 to 20, because 
the nonmagnetic GS is massive. In LL2, however, the 
correct magnetic field H for to cannot be obtained by 
j because the excitation with SM ± 1 has a finite gap 
proportional D. Thus we have to use another form 

H$ = [E(L, M + 2) - E(L, M - 2)]/4 = H + o(i), (12) 

where only even M is available. H c \ and H C 2 at the 
boundary of LL2 are estimated by the Shanks transfor- 
mation applied to (E(L, 2) - E(L, 0))/2 and (E(L,L) - 
E(L — 2))/2, respectively. Fig. ^| shows the completed 
H-D phase diagram. 
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FIG. 5. H-D phase diagram. H C 3 (a dot-dashed line) and 
H C 4 (a thick solid line) were obtained by plots of (D c , H^) 
for various m with L = 20 and 18, and a curve fitting. H c i (a 
thick dashed line) and H C 2 (a long dashed line) were obtained 
by the Shanks transformation applied to Ai (LL1) and A2 
(LL2) corresponding to m — 0+ and m = 1— , respectively. 
The field induced transition from LL2 to LL1 occurs at H C 3 
fon —1.49 < D < —0.35, and the re-entrant transition also 
exists for -1.49 < D < -1.32. 



VI. MAGNETIZATION CURVE 

Finally we investigate the behavior of the magnetiza- 
tion curve, in particular around the critical fields. If we 
assume the form around the critical field H c correspond- 
ing to the magnetization m c as 

m-mc-iH-Hc) 1 / 5 , (13) 

the exponent 5 can be estimated by the form jllj 

'"(j^l/M^W (1 ^>' (14) 

where f(L) is given by 

f(L) = E(L, M c + 2) + E(L, M c ) - 2E(L, M c + 1), (15) 

where M c — m c L. For LL2 the forms need some mod- 
ifications like L + 2 -» L + 4, M c + 2 -» M c + 4 and 
M c + 1 M c + 2. The method applied to H cl at L> = 
(isotropic case) gave the result 5 = 1.9 ± 0.1 consistent 
with 5 = 2 predicted by some effective Hamiltonian the- 
ories. ]ll[ Our present application of the method to H c i 
at D = -2.0 in LL2 yields 5 = 2.03 ± 0.04 which sup- 
ports the equivalence between the present model (Q) and 
the anisotropic S = 5 chain, discussed above. To esti- 
mate 5 around H c3 , we applied it to the phase bound- 
ary D c = -1.44 for m = ±, and got 6 = 0.9 ± 0.3 
which suggests 5=1. It implies that the magnetization 
curve is linear without any discontinuity around H c ^ , be- 
cause any jump at H c ^ whoud make the estimated 5 very 
small. However, the possibility of a change in the deriva- 
tive of the magnetization curve at H C 3 still exists. To 
test it, we give the ground state magnetization curve for 



D = -1.0, -1.44 and -2.0 in Fig. g, using the form 
(|ll]) in LL1 and (|l|) in LL2 for L = 20. The curve for 
D = -1.0 has H c3 = -1.59, and the one for D = -1.44 
has H c3 = -2.27 and H c4 = -2.58. In Fig. | solid 
and dashed lines stand for LL1 and LL2, respectively. 
At least within this analysis, no significant anomalous 
behavior can be detected around H c3 or H c4 . 
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FIG. 6. Magnetization curve in the ground state for 
D = —1.0, —1.44 and —2.0. Solid and dashed lines corre- 
spond to LL1 and LL2, respectively. The transition from 
LL2 to LL1 appears at H C 3 = —1.59 for D = —1.0, and at 
H C 3 = —2.27 for D = —1.44. The re-entrant transition oc- 
curs at // C 4 = —2.58 for D = —1.44. No significant anomalous 
behavior appears at any critical fields. 

VII. DISCUSSIONS 

To consider the relation of the transition at H c3 to the 
spin flopping, the line of the spin flopping in the classical 
limit with (Sj) 2 = 1, which is D = -2m 2 /(l + m 2 ), is 
shown together in Fig. |4j In the classical spin flopping m 
jumps from to the line and the magnetization process 
is equivalent to the Ising model for D < — 1. In a way 
the transition at H c3 is similar to the spin flopping be- 
cause the transverse spin correlation has an exponential 
decay for LL2, while power-law one for LL1, which im- 
plies that even infinitesimal interchain interaction would 
yield the canted Neel order in LL1. LL2, however, is es- 
sentially produced by quantum effect, with no classical 
correspondence. 

The quasi-lD antiferromagnet CsNiChj can be a good 
candidate to test the transition at H c3 because it is in 
the Neel phase in GS due to interchain interaction, p^ ] 
Thus it is interesting to consider the quasi-lD case. The 
interchain interaction enhances the most dominant spin 
correlation and leads to the long-range order. Thus it 
is expected that LL1 is changed into the canted Neel 
order, while LL2 into the incommensurate SDW with k = 
(1 — m)n along H . In addition the transition at H c3 could 
become first-order one with a jump in the magnetization 
curve. 
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VIII. SUMMARY 



The finite cluster calculation and some size scaling 
analyses suggested that a field induced phase transition 
between two different Luttingcr liquids occurs at H c z , be- 
tween H c i and H c2 , in the ID S = 1 antifcrromagnet with 
easy-axis anisotropy for —1.49 < D < —0.35. In addition 
a re-entrant transition can appear for —1.49 <D< —1.32. 
The transition is reminiscent to the spin flopping, but 
it gives rise to no significant anomalous behavior in the 
magnetization curve. 
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